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Abstract. We briefly review the properties and physical consequences of quasiperiodic 
oscillations (QPOs) seen on many occasions in the X-ray emission from black-hole binary 
systems. High frequency QPOs [y > 40 Hz) continue to be scrutinized as effects of general 
relativity, with new attention to the role of resonances in their formation. Low-frequency 
QPOs (0.05 to 30 Hz) exhibit complicated behavior, with occasions of high amplitude 
and particular correlations with some X-ray spectral parameters. QPO mechanisms are 
a requirement for any physical model seeking to explain either (1) the non-thermal X-ray 
spectrum that is commonly seen and is usually stronger than the accretion disk at times 
of highest luminosity, or (2) the hard X-ray spectrum evident when there is a steady type 
of radio jet. 



1. Introduction 

Most of the brightest celestial X-ray sources recorded in a given year are transient 
outbursts from black hole binaries in the Galaxy. The eruptions are understood 
as a consequence of a low rate of mass accretion from the companion star [1]. 
The material gradually fills the outer regions of an accretion disk surrounding 
the black hole. When the disk surface density reaches a critical value, matter 
spirals into the inner disk where it reaches X-ray emitting temperatures before 
falling into the black hole event horizon [2]. Such X-ray novae are the parent 
population for identifying black holes that are remnants of massive stars. The 
mass of the compact object (typically 5-15 M & ) is deduced from radial velocity 
studies of the companion star. In most cases such measurements are only possible 
when the system has returned to relative quiescence, and the companion can be 
seen against the glare of the hot gases in the disk. In establishing the nature of 
the compact object, the critical argument is whether the mass exceeds the upper 
limit (~ 3.0 Mq ) for the mass of a neutron star. There are now 17 "dynamical 
black hole" binaries in the Milky Way or LMC: 14 were first noticed as bright 
X-ray novae [3], while the other 3 cases are persistent X-ray sources with O/B 
type companions. 

Observations with the Rossi X-ray Timing Explorer (RXTE) have pioneered 
efforts to further study black holes and their occasional relativistic jets via broad- 
band X-ray observations during active states of accretion. The X-ray timing and 
spectral properties convey information about physical processes that occur near 
the black hole event horizon, and one of the primary research goals is to obtain 
constraints on the black hole mass and spin using predictions of general relativity 
(GR) in the strong-field regime. There is also the need to understand accretion 
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physics for each of the four distinct emission states that are displayed by so many 
accreting black holes systems. In this paper we describe recent advances in these 
topics, with particular attention to the diverse forms of quasiperiodic oscillations 
(QPOs) [4] that have been detected during the extensive monitoring programs 
for X-ray transients conducted with RXTE. 

2. High Frequency Oscillations in Accreting Black Hole Systems 

The topic of high-frequency QPOs (HFQPOs) in black hole binaries (40-450 Hz) 
has evolved substantially in the last year. Transient HFQPOs are detected with 
RXTE in 5 sources (4 dynamical black holes and 1 candidate) . These subtle oscil- 
lations have rms amplitudes that are typically only ~ 1 % of the mean count rate. 
The associated power spectra are shown in Fig. [l], and their properties are summa- 
rized in Table 1. Three of these sources exhibit pairs of HFQPOs: GRO J1655-40 
(300, 450 Hz; [5]), XTE J1550-564 (184, 276 Hz; [6, 7]), and GRS 1915+105 (40, 
67; [8]). Furthermore, in Section 3 below we report the detection of additional 
HFQPOs in GRS 1915+105 at 164 and 328 Hz. 

The commensurate frequencies (3:2 ratio) in pairs of HFQPOs from 
GRO J1655-40 and XTE J1550-564 [9] may suggest a cause in some type of 
resonance phenomenon involving oscillations proscribed by general relativity, as 
originally proposed by Abramowicz & Kluzniak [10]. The observations cannot 
be interpreted as non-sinusoidal structures in the wave pattern, in the sense of 
harmonics in a Fourier series, because the individual detections (in a given energy 
band) appear as a single peak in the power spectrum, and their frequencies are 
recognized as 2v or 2>v only when the ensemble of results is examined. 



TABLE 1. Table 1: High Frequency QPOs in Black Hole Systems 



Source # RXTE QPO Freq. Coherence Detections Energy Ampl. 



Name 


Obs. 


Hz (±) 


(v/FWHM) 


# (3 a) 


(keV) 


rms % 


GRO J1655-40 


81 


300 (23) 
450 (20) 


4.4-11.6 
8.8-14.8 


7 
6 


2-30 
13-30 


0.5-1.0 
2.5-5.3 


XTE J 1550-564 


364 


184 (26) 


2.6-9.1 


10 


2-30 


0.6-1.7 






272 (20) 


4.7-14.3 


14 


6-30 


3.6-6.9 


GRS 1915+105 


448 


67 (5) 
41 (1) 
164 (2) 
328 (4) 


3.2-22.0 
5-20 
5-7 
14-16 


29 
5 

2 sums[14] 
2 sums [14] 


2-30 
13-30 
13-30 
13-30 


0.4-1.8 
1.8-3.2 
2.0-2.2 
1.0-1.3 


4U 1630-47 


301 


184 (5) 


5-9 


sum [30] 


6-30 


1.0-1.2 


XTE J1859+226 


135 


193 (4) 


3-5 


sum [48] 


6-30 


1.7-2.0 



The resonance hypothesis [10] has been discussed in terms of accretion blobs 
following perturbed orbits in the inner accretion disk. Unlike Newtonian gravity, 
GR predicts independent oscillation frequencies for each spatial coordinate for 
orbits around a rotating compact object, as seen from the rest frame of a distant 
observer. GR coordinate frequencies and their differences (i.e. beat frequencies) 
have been proposed to explain some of the X-ray QPOs from both neutron stars 
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Figure 1. HFQPOs in five black hole binary systems. The energy band is 6-30 keV unless 
otherwise indicated. These subtle oscillations are only visible during a fraction of the observations 
for each source, and details are given in Table 1. 



and black holes [11]. At the radii in the accretion disk where X-rays originate, the 
GR coordinate frequencies are predicted to have varying, non-integral ratios. The 
commensurate frequencies in pairs of HFQPOs, noted above, can therefore be in- 
terpreted as a signature of resonance, e.g. between pairs of coordinate frequencies 
in the inner disk. Unlike the azimuthal and polar coordinates frequencies, the ra- 
dial coordinate frequency reaches a maximum value and then falls to zero as the 
radius decreases toward the location of the innermost stable orbit (see [12], [13], 
and references therein). This ensures the possibility of commensurate coordinate 
frequencies somewhere in the inner disk. For example, there is a wide range in 
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the dimensionless spin parameter, a*, where one can find a particular radius that 
corresponds to a 2:1, 3:1, or 3:2 ratio in the orbital and radial coordinate fre- 
quencies. A resonance between the polar and radial coordinate frequencies is also 
possible. In the resonance concept, nonlinear perturbations may grow at these 
radii, ultimately producing X-ray oscillations that represent some combination 
of the individual resonance frequencies, their sum, or their difference. However, 
there remain serious uncertainties as to whether such structures could overcome 
the severe damping forces and emit X-rays with sufficient amplitude and coher- 
ence to produce the observed HFQPOs [14]. 

Models for "diskoseismic" oscillations adopt a more global view of the inner 
disk as a GR resonance cavity [15, 16]. This paradigm has certain attractions for 
explaining HFQPOs, but integral harmonics are not predicted for the three types 
of diskoseismic modes derived for adiabatic perturbations in a thin accretion disk. 
Clearly, there is also a need to investigate the possibility of resonances within the 
paradigm of diskoseismology. 

As initially noted by Strohmayer [5], the HFQPO measurements may be seen 
to suggest substantial black hole spin. For the 3 sources with pairs QPOs, op- 
tical or IR studies have yielded measurements of the black hole mass. In each 
case the fastest QPOs exceed the maximum orbital rotation frequency {y^) at 
the innermost stable circular orbit around a Schwarzschild black hole [17] (i.e. 
dimensionless spin parameter, a* = 0). If the maximum Keplerian frequency is 
the highest frequency at which a QPO can be seen, then the measurements re- 
quire prograde spin, e.g. a* > 0.15 for GRO J1655-40. On the other hand, even 
higher values of the spin parameter are required if the QPOs represent either a 
resonance in coordinate frequencies [7] (a* > 0.25), or pairs of diskoseismic modes 
that exhibit commensurate frequencies by chance [18] (a* > 0.9). 

Encouragement for linking HFQPOs and GR oscillations can be found in the 
fact [9] that the HFQPOs in XTE J1550-564 and GRO J1655-40 scale ~ M -1 . 
This result is generally consistent with the known mechanisms related to disk 
oscillations in the strong-field regime of general relativity, as long as the values 
of the spin parameter (a*) are similar for these two black holes. These results 
illustrate both the quantitative value of HFQPO detections as a means of probing 
the physical properties of black holes, and also the need to continue efforts to 
independently measure black hole masses via dynamical optical studies. The spin 
measurements would be of further value towards evaluating the role of black hole 
rotation in the production of jets [19]. 

Finally, the investigation of HFQPOs and the the energy spectra for both 
XTE J1550-564 and GRO J1655-40 show identical patterns of behavior that 
demand explanation [9]. There is a systematic increase in the strength of the 
hard X-ray power-law component (i.e. the type with index ~ 2.5) as the observed 
HFQPO shifts from 3v to 2v. We also note that HFQPOs are generally not 
detected either when the accretion disk dominates the spectrum (i.e. power-law 
is absent), nor when there is radio emission and the power-law is flatter (index 
~ 1.7). These patterns of spectral behavior were used to select the data that yields 
a new HFQPO detection for GRS 1915+105, as described in the next section. 

3. New HFQPOs with Harmonics in GRS 1915+105 

GRS 1915+105 is the most extraordinarily variable X-ray binary known [20]. It is 
also the first member of the "microquasars" , a handful of Galactic X-ray sources 
that produce relativistic jets [21], as seen in the "super luminal" separation of 
bipolar knots in sequences of radio images. The system is a confirmed black hole 



Workshop Proceedings 



5 



binary that contains a ~ 14 M black hole [22]. A rich variety of correlations 
have been found between X-ray, radio, and infrared outbursts [21, 23, 24, 25, 26, 
27], contributing substantially to the overall success in linking mass ejections to 
conditions in the inner accretion disk. 

There is an extensive archive of ~ 800 pointed observations of GRS 1915+105 
obtained with RXTE during the last six years. About half of these show sub- 
stantial variability, while the others show roughly steady flux with a/fj, < 0.18, 
evaluated in 1 s bins. Studies of the energy division between the thermal/disk 
and hard power-law components suggest that most of the steady states are either 
hard states with persistent radio emission or soft states in which most of the flux 
originates in the disk [28]. The transient HFQPO at 67 Hz [29] tends to appear 
in the latter condition, especially when the spectrum is soft and very bright [30]. 
On the other hand, the steady states of GRS 1915+105 offer limited exposures 
to conditions when both the disk and power-law exhibit high luminosity, which 
corresponds to the times when harmonic HFQPOs appear in XTE J 1550-564 
and GRO J1655-40. We are therefore motivated to investigate the variable states 
of GRS 1915+105, using a scheme to sort the exposures into quadrants of the 
color-intensity diagram, in order to accumulate significant exposure time when 
the source is bright and moderately hard. 



> 



X 

< 
u 

PL 



0.08 



0.06 - 



0.04 



0.02 



GRS1915+105 Sept. 1997 

H3 




S3 



2000 4000 6000 8000 
Count Rate (c/s/PCU) 



I 4 



Figure 2. Color-intensity diagram for GRS 1915+105 using RXTE PCA observations on 10 
days during 1997 September 5-29. The data from all operating PCUs have been normalized and 
combined For each 16 s exposure. The dashed lines illustrate the sorting boundaries used to 
extract times when the source has different levels of the intensity and the relative strength of 
the power-law component. 

For the present analysis, we select the 14 RXTE observations during 1997 
September 5-29. These data offer long exposures, some rare bright-hard-steady 
states, and types of wild variations that include frequent excursions into bright- 
hard states. The latter cases are highlighted by light curves of the 9 classification 
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Figure 3. Average power spectrum for groups H2 and H3 (see Fig. 2), selecting only the PCA 
detections above 13 keV. The arrows indicate weak features at 164 Hz and the locations of the 
next harmonics (i.e. 328 and 492 Hz). 

of Belloni et al. [20]. The background-subtracted PCA spectra are analyzed at a 
time resolutions of 16 s. The output of all operating PCUs is combined for each 
time interval, after normalizing each PCU to the characteristics of PCU #2. The 
PCU normalization is defined using six energy intervals and simply scaling the 
mean count rate per interval to the results for PCU #2, using the mean spectrum 
derived for the subset of the observations in which both PCUs were operating. 
Locations on the color-intensity diagram are computed using a hardness ratio 
(HR) of normalized source counts at 13-30 keV vs. 2-13 keV. The resulting color- 
intensity diagram and the data sorting scheme are shown in Fig ||. The hard 
("H") and soft ("S") regions are separated by the line, HR = 0.045, and there 
are three intensity levels (1,2,3), as shown in the Figure. 

The average power density spectra (PDS) are then computed for the data in 
each of the six regions of the color- intensity diagram. The search for HFQPOs 
(see methods in [9]) produces detections restricted to regions H2 and H3, using 
data confined to the energy range of 13-30 keV. These PDS are shown in Fig |f 
The H2 PDS shows a peak with a significance of 5.7 a at 163 ± 5 Hz with an rms 
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Figure 4- Average power spectrum for groups H2 and H3 (13-30 keV) with the best fit for QPO 
features and the power continuum. The QPO at 164.9 Hz has a significance of 9 a. The QPO 
at the first harmonic (329.8 Hz) is detected at 3.7 a, while the hint of a second harmonic (484.8 
Hz) is not statistically significant (1.5 a). 



amplitude of 1.56 ± 0.14 % of the mean count rate at 13-30 keV. Similarly, the 
H3 PDS yields a 4.4 a HFQPO at 166 ± 3 Hz with amplitude 1.73 ± 0.20 %. 

We then combine the H2 and H3 PDS (weighted mean values per bin) and 
fit the results for an HFQPO with harmonics. The data and best-fit model are 
shown in Fig ||, and the HFQPO parameters are included in Table 1. With the 
first harmonic included in the model, the HFQPO at 165 ±3 Hz has a significance 
above 9 a. The second peak at 330 Hz has a significance of 3.7 a, while there is 
a hint (not significant) of the next harmonic at 495 Hz. 

In conclusion, the spectral characteristics of HFQPOs in XTE J1550-564 and 
GRO J1655-40 have provided guidance for finding new HFPQOs at 165 and 330 
Hz in GRS 1915+105. The 165 Hz QPO in GRS 1915+105 then appears closely 
related to the QPO systems with v = 150 Hz and 92 Hz in GRO J1655-40 
and XTE J1550-564, respectively. Given the optical mass estimates for all three 
black holes, the results suggest faster spin in GRS 1915+105, compared to the 
other two sources, if the QPOs originate from a common inner disk oscillation 
related to general relativity. On the other hand, the QPOs at 41 and 67 Hz in 
GRS 1915+105 may arise from a different physical mechanism. 

4. Low Frequency QPOs and Problems in Accretion Physics 

The X-ray PDS of many black hole transients display low frequency QPOs (LFQ- 
POs), roughly in the range of 0.05 to 30 Hz. The oscillations can be particularly 
strong, with peak to trough ratios as high as 1.5 [31]. LFQPOs have been traced 
out to energy bands as high as 60-124 keV [32]. The measurement properties of 
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LFQPOs are summarized in Table 2. Some authors pursue separate interpreta- 
tions for HFQPOs and the more common LFQPOs because they seem to evolve 
separately and may sometimes coexist [7]. However, others argue that there are 
unification schemes involving both neutron stars and black holes in which partic- 
ular mechanisms may operate across a wide range in frequency [33]. 

Investigations of LFQPO properties versus spectral parameters of the disk 
and power law components generally show that the LFQPO frequency is most 
often correlated with the disk flux [34, 35]. However, the rms amplitude generally 
increases between 2 and 20 keV, and the amplitude spectrum therefore resembles 
the power-law component rather than the accretion disk [32]. This dualism in 
correlations encourages LFQPO investigations motivated to deduce the origin 
of the power-law spectrum. One major complication here is that broad-band 
spectral studies tend to emphasize the importance of distinguishing two forms 
of X-ray power law spectra [36]. At highest luminosity (the "very high" state), 
the spectrum is dominated by a power law with a photon index, T ~ 2.5, which 
may extend to 500 keV or higher [37]. Alternatively, in the "low/hard" state, the 
power-law component shows a flatter spectrum (T ~ 1.7) that is cut off above 
~ 100 keV, and this state is almost always associated with a steady radio jet [38]. 
On the other hand, the PDS and QPO properties of these two states may appear 
rather similar [39, 31]. 



TABLE 2. Table 2: Properties of Low Frequency QPOs in Black Hole Systems 



Property 


Value 


Comments 


Frequency range 


0.05 - 30 Hz 


most 1-7 Hz 


rms amplitude 


1 - 20 % 


2-30 keV; higher at 6-30 keV 


Q = v/Av F whm 


3-30 


typical ~ 8.5 


Phase Lag 


-0.1 - 0.2 


2-6 keV vs. 13-30 keV 


Harmonics 


sometimes 2 and 0.5 ^ 


e.g. XTE J1550 and GRS 1915 


In Low-Hard State 


often 


Y pl ~ 1.7 + radio flux 


In High-Soft State 


no 


thermal spectrum, weak power law 


In Very High State 


yes 


defining item; Ypl ~ 2.5 


Min. Turn-On Timescale 


seconds 


GRS 1915+105 light curves 


Max. Longevity Timescale 


months 


GRS 1915+105 low-hard state 


Freq. Correlation 


with Disk Flux 


all other params worse 


Ampl. Correlation 


with Photon Energy 


like power-law spectrum, not disk 



For the strongest QPOs in GRS 1915+105, efforts were made to track the 
individual oscillations to determine the origin of the frequency drifts and to mea- 
sure the average "QPO-folded" oscillation profile [29]. The results show a random 
walk in QPO phase, with only minor departures from a sinusoidal waveform. The 
ramifications of this study for QPO models remain uncertain. There are now a 
large number of proposed LFQPO mechanisms in the literature, and a full review 
is beyond the scope of this paper. The models include global disk oscillations [40] 
as well as oscillations in the radial position of complex accretion structures, such 
as shock fronts [41] or a transition layer between the disk and a hotter Comp- 
tonizing region [42]. The construction of complex accretion models is driven not 
only by the need to account for the power-law spectrum, but also by the relatively 
slow frequencies of LFQPOs. Given the capacity of LFQPOs for high amplitudes 
and long-term stability (see Table 2), it would be attractive to associate them 
with Keplerian frequencies in the disk; however, the relevant disk radii are well 
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beyond the expected limits of the X-ray emitting region. An alternative model 
known as the "accretion-ejection instability" invokes spiral waves in a magnetized 
disk [43] , with a transfer of energy out to the radius of corotation with the spiral 
wave. This model uses magnetic field effects to account for QPOs, the power-law 
spectrum, and perhaps relativistic outflows as well. Further details are given in 
other papers of this workshop. 

Recent efforts to understand LFQPOs have given attention to the phase lags 
associated with these oscillations and their harmonics. The analysis technique 
uses Fourier cross spectra to measure both the phase lags and the coherence 
parameter (versus frequency) between different X-ray energy bands, e.g. 2-6 
vs. 13-30 keV. Unexpectedly, both positive and negative phase lags were found 
[44, 45, 46, 39], and suggestions were made to classify LFQPOs by phase lag 
properties [44]. While the expansion of LFQPO subtypes may not be a welcome 
complication, it has been shown that recognition of the phase lag properties al- 
lows us to better understand the relationship between certain types of QPOs 
and the conditions in the disk [7]. This is illustrated in Fig 0. Furthermore, data 
selections by LFQPO subtypes provided guidance that helped to establish the 
harmonic nature of HFQPOs in XTE J1550-564 and the relationship between 
LFQPOs and HFPQOs in that source [9]. 



XTE J1550-564 1998-1999: LFQPO Freq. vs. Disk Flux 
12 




0.1 0.2 0.3 

Disk Flux (2-20 keV) 



Figure 5. Correlation between LFQPO frequency and the apparent disk flux for XTE J1550-564 
during its 1998-1999 outburst. The plotting symbols distinguish the QPO type: Type A (broad 
QPOs with phase lags in soft X-rays) - open triangles, Type B (narrow QPOs with hard lags) 

- open squares, Type C (small lags and strong harmonics) - filled circles, and anomalous QPOs 

- 'x'. The correlation between these quantities is only demonstrated for the C type LFQPOS. 

In closing, it must be noted that the relevance of X-ray timing studies requires 
discussions of broad power peaks in the PDS [47, 48], which are features wider 
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than the limit (Q ~ 2) used to define QPOs. Broad power peaks are involved 
in important studies such as the means to distinguish accreting black holes from 
neutron stars [49]. The evolution of broad power peaks has also been linked to 
major behavioral changes in a black hole binary, as shown in disappearance of 
a broad power feature just at the time when Cyg X-l recently lost its ability to 
remain in the low-hard state and maintain a steady jet [51]. 
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